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I. INTRODUCTION 

Processes with identified final-state hadrons play im- 
portant roles in QCD. Foremost, they provide crucial 
information on fragmentation functions and hence, ul- 
timately, the hadronization mechanism. Modern anal- 
yses [IH1] of fragmentation functions variously use data 
for single-inclusive annihilation (SIA) e + e~ — > hX, semi- 
inclusive deep-inelastic scattering (SIDIS), lp — > £hX, 
and pp — > hX, where h denotes a final-state hadron. 
Hadron production observables also serve as powerful 
probes of nucleon or nuclear structure. In particular, 
SIDIS measurements with polarized beams and/or tar- 
gets have by now become indispensable tools for inves- 
tigations of the spin structure of the nucleon in terms 
of helicity parton distribution functions and transverse- 
momentum dependent distributions [5J. Finally, hadron 
production data also test some of our key concepts in the 
theoretical analysis of QCD at high energies, among them 
factorization, universality, and perturbative calculations. 

In the present paper, we address higher-order pertur- 
bative corrections to two of the key hadron production 
processes, SIA e+e~ -> hX and SIDIS lp -> IhX. Our 
study is very much motivated by the recent advent of 
data for these reactions with unprecedented high preci- 
sion. The BELLE collaboration at KEK has presented 
preliminary data [6] for pion and kaon multiplicities in 
SIA that extend over a wide range in values of the frag- 
mentation variable xe — 2E h /y/s, where E h is the en- 
ergy of the produced hadron in the e + e" center-of-mass 
system, and y^i = 10.52 GeV the collision energy. The 
BELLE data cover the region 0.2 < xe < 0.97, with a 
very fine binning and extremely high precision partly at 
the sub-1% level. New preliminary high-statistics SIDIS 
data have been shown by the HERMES \7[ and COM- 
PASS j8] collaborations over the past year or so. 

In the kinematic regimes accessed by these experi- 
ments, perturbative-QCD corrections are expected to be 
fairly significant. In case of e + e~ — > hX at BELLE, as 
xe increases toward unity, the phase space for real-gluon 
radiation is very restricted, since most of the initial lep- 
tonic energy is used to produce the observed hadron and a 
recoiling unobserved final state. When this happens, the 
infrared cancellations between virtual and real-emission 



diagrams leave behind large logarithmic higher-order cor- 
rections to the basic e + e~ — > qq cross section. These 
logarithms are very similar in nature to the "threshold 
logarithms" encountered in hadronic scattering processes 
when the energy of initial partons is just large enough 
to produce a given final state. Near xe = 1, it then be- 
comes necessary to take the large corrections into account 
to all orders in the strong coupling, a technique known 
as threshold resummation. The SIDIS measurements, on 
the other hand, are characterized by two scaling vari- 
ables, Bjorken-x and a variable z given by the energy of 
the produced hadron over the energy of the virtual pho- 
ton in the target rest frame. The cross section is typically 
defined to be differential in both. Large logarithmic cor- 
rections to the SIDIS cross section arise when either one 
of the corresponding partonic variables becomes large. 
The most important effects arise when both are large, 
which is typically the case for the presently relevant fixed- 
target kinematics. As we shall discuss, in this case the 
logarithmic terms can be simultaneously resummed to all 
orders within the threshold resummation framework. 

Previous work [5] has established a close correspon- 
dence between threshold resummation for the Drell-Yan 
process, double-inclusive annihilation e + e~ — > hih 2 X, 
and a variant of SIDIS for which one considers the cross 
section differential in the product xz of the two scaling 
variables mentioned above, rather than in each of them 
separately. For this set of observables, the structure of 
the threshold logarithms turns out to be identical, up to 
trivial differences in the hard scattering functions that 
multiply the logarithms. One therefore can derive the re- 
summation for each of the processes in the same manner, 
using exponentiation of eikonal diagrams in color-singlet 
processes [101 lllj . This approach was termed crossed 
threshold resummation in [9] and will be the framework 
for our present analysis of SIDIS. 

Resummation for e + e~ — > hX was addressed in detail 
in Refs. [T2TU5| . In [T2] the next-to-leading logarithm 
(NLL) expressions were presented, while extensions to 
next-to-next-to-leading logarithm and even next-to-next- 
to-next-to-lcading logarithm were provided in Refs. [T3] 
and [2], respectively. The global analysis of fragmen- 
tation functions of Ref. [3] in fact includes NLL thresh- 
old resummation effects for e + e~ —¥ hX for the lower- 
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energy data. They are found to improve the theoretical 
description and the quality of the fit to the data. Given 
the very high precision of the new BELLE data, a new 
phenomenological analysis of resummation effects is now 
timely and will be presented in this work. For now, we 
will restrict ourselves to resummation at NLL, which cap- 
tures the main effects. As for SIDIS, an expression for 
the NLL resummed cross section was stated in [T2] , which 
turns out to be analogous to the related expressions 
for the rapidity-differential Drell-Yan cross section in 
terms of double-Mellin moments (see also [MHU]). Phe- 
nomenological studies of resummation effects in SIDIS 
have never really been presented in the literature, ex- 
cept briefly in 18J. In the present paper, we derive the 
NLL resummed expression for SIDIS, making use of the 
techniques of [9]. We also present numerical results as 
relevant for comparisons to the recent SIDIS data. 

Sections HI1 and Mil collect our technical derivations. 
Since the resummation for SIDIS is the main new re- 
sult, we mostly focus on this process and only briefly 
review the well-known results for e + e _ — > hX and inclu- 



SIDIS hadron multiplicities are defined by 



sive deep- inelastic scattering (DIS). Section IV presents 
our phenomenological results. 



SIDIS 



d 3 cr h /dxdydz 
d 2 a /dxdy 



(3) 



where d 2 a/dxdy is the cross section for inclusive DIS, 
Ip — > IX, given by 
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dxdy Q 2 



1 + (1 - y) s 



F T {x,Q 2 )+ ] -^F L (x 1 Q 2 ) 

y 

(4) 

Here Ft = 2Pi and Fl = F^/x, with the standard inclu- 
sive structure functions F\ , Fl ■ Usually, the numerator 
and denominator of ^ are averaged over suitable bins 
in x and y. In order to investigate higher-order effects 
on SIDIS multiplicities, we have to consider QCD cor- 
rections to both the SIDIS and the inclusive DIS cross 
section. Since the latter has been treated very exten- 
sively in the literature, we will focus here on d 3 a h /dxdydz 
and only briefly summarize some of the known results for 
d 2 a I dxdy. 
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A. SIDIS cross section at next-to-leading order, 
and Mellin moments 



II. RESUMMATION FOR SIDIS 
MULTIPLICITIES 

We consider semi-inclusive deep-inelastic scattering, 
l(k)p(P) — > £(k')h(Pf l )X , where we have indicated the 
momenta of the involved particles. The momentum q of 
the highly virtual photon exchanged between the incom- 
ing electron and proton is given by q = k — k'. We define 
the usual variables 

Q 2 = -q 2 = -(k-k') 2 , 

Q 2 
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We have Q 2 — xys, with yfs the ip center-of-mass energy. 
In the current fragmentation region that we will consider 
here, the SIDIS cross section may be written as [I9H24] : 
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Using factorization, the transverse and longitudinal 
structure functions F% = 2F/ 1 and T\ = F£/x in Q 
are given by (i = T, L) 
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1 dx f 1 dz 



x C f'f 



X 7 z, 
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where /(£,M 2 ) denotes the distribution of parton / = 
q, q, g in the nucleon at momentum fraction £ and scale 
/i, while Dj, (£, /x 2 ) is the corresponding fragmentation 
function for parton /' going to the observed hadron h. 
For simplicity, we have set all factorization and renor- 
malization scales equal and collectively denoted them by 
fx. The hard-scattering coefficient functions Cj-,^ can be 
computed in perturbation theory: 



L f'f ~ U f'f 



a s (y 2 ) 
2tt 



0(a 
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where, again, i = T,L. To lowest order (LO), only the 
process j*q —¥ q contributes, and we have 



C^°\x,z) = e 2 q S(l-x)S(l-z), 
C^°\x,z) = 0, 



(7) 



where a is the fine structure constant. F? and F£ are the 
transverse and the longitudinal structure functions; they 
are related to the more customary structure functions F^ 
by T% = 2Ff and F\ = F£/x. 



and F£ 



with the quark's fractional charge e q . Beyond LO, also 
gluons in the initial or final state contribute. The full set 
of the first-order coefficient functions 
collected in the Appendix. 



Cf'f are 
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Since threshold resummation can be derived in Mellin- 
moment space, it is useful to take Mellin moments of 
the structure functions IFj, and J 7 ^. Since x and z are 
independent variables, we take moments separately in 
both [HdS]. We define 



P*{N,M,Q 2 )= ( dxx N - 1 f dzzJ 1 - 1 ^(x,z,Q' 
Jo Jo 



We then readily find from (JsJ 

/,/' 
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(9) 



where 



f N (n 2 ) = / dxx^fix,^), 
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dxx 



C'f>f (x,z, ^,a s (fi 2 ) 



(10) 

Thus, the Mellin moments of the structure functions 
are obtained from ordinary products of the moments 
of the parton distribution functions and fragmentation 
functions, and double-Mellin moments of the partonic 
hard-scattering functions. For the perturbative expan- 
sion given in J6J), we have for the latter in lowest oder 
according to (Fff: 



C^°\N,M) = e\, 
C^°\N,M) = 0. 



(11) 



The corresponding moments of the next-to-leading or- 
der (NLO) terms [25] are also provided in the Ap- 
pendix. 



Choosing for simplicity fip — Q, we have 



Cj-W(:M) ~ e 2 q C F 



(1 - x)5(l - z) 



+25(1 -£)[—, r 1 ) +25(1-2)' 



V i - z 



1-x 



(l-x) + (l-z) + 
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corresponding in moment space to 



CT>W(N,M) 



elC F 



TV- 



(lniV + lnM)' 



(13) 



where N = Ng^ b , M = Me 7B , with ^ E the Euler con- 
stant. Here we have only kept contributions that are nei- 
ther suppressed as 1/N, nor as 1/M in moment space. 
The terms given in ( 12 ) therefore always contain two dis- 
tributions, one in x and one in z. 

Threshold resummation addresses the logarithms in 
N and M to all orders in the strong coupling con- 
stant a s . More precisely, it captures terms of the form 
a k s In™ N\n m M, with n + m < 2k. We now discuss the 
derivation of the resummed expression for the SIDIS co- 
efficient function C qq A >(N,M). Since the leading-order 
process is 7*g — > q scattering, and since both the incom- 
ing and the outgoing quark are "observed" , the treatment 
has much in common with that for the total Drell-Yan 
cross section, or for its "crossed" versions considered in 
Ref. [S]. A significant difference is, however, that in the 
present case two independent Mellin moments, N and 
M, have to be considered. At large N and M, or equiva- 
lently x and z, all gluon radiation from the basic process 
j*q — > q becomes soft, since we have the relation [§] 



(l-x) + (l-z) 



2k° 



(14) 



B. Resummation of the SIDIS coefficient function 



As one can see from Eq. ( 49 ) , the NLO coefficient func- 
tion Cqq (x, z) receives large corrections near x, z 1. 



where k° is the total energy of gluon radiation. The 
coefficient function may then be evaluated in the eikonal 
approximation for the quarks and/or antiquarks involved 
in the hard scattering. In moment space, the eikonal hard 
scattering functions exponentiate, leading to [9"HTT] 
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(15) 



which is valid to next-to-leading logarithmic (NLL) ac- 
curacy. Here, A q (a s ) is a perturbative function: 



A q (a s ) = °±AU + 

7T \ 7T 



with 



A (2) 



67 _ 7T^ 
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(16) 
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(17) 

where Cf — 4/3, Ca = 3 and iV/ is the number of active 
flavors. 

Up to corrections that are exponentially suppressed at 



large N, M, the integral over £ in ( 15 1 can be carried out 
analytically, and one finds 



C iQ 2 



IniV + lnM 



hi 



Q 



, (18) 



where Kq is a Bessel function. It arises when we extend 
the £ integral to < £ < oo. It is instructive to confront 



this with the analogous expression for the Mellin-iV mo- 
ments of the partonic total Drell-Yan cross section, which 
reads El [11]: 



1 d A 
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From this comparison one can immediately see that the 
result for the resummed SIDIS cross section can be ob- 
tained from the one for the total Drell-Yan cross section 
by simply setting N — > VNM. In the case of SIDIS, the 
moments N and M independently fix the light-cone plus 
component £ and the minus component £ = fc^/(£Q 2 ) of 
the soft gluon momentum, resulti ng in the slightly more 
elaborate form of the exponent in ( 18 ). Likewise, the MS- 



subtraction of collinear divergencies in the initial and the 
final state in SIDIS yields the terms IniV and InM, re- 
spectively whereas in Drell-Yan one has the contribution 
2 In ./V from the two initial partons. 

With (18), the final resummed coefficient function be- 



comes in the MS scheme: 
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(20) 



Here we have included a perturbative function H qq that 
collects the hard virtual corrections to 7*5 — > q scatter- 
ing. For resummation at NLL, one needs to know H qq to 
first order in the strong coupling, which may be derived 
by expanding (20) to 0(a s ) (keeping only logarithmic 
terms in the exponent) and comparing to the explicit 
NLO expression (13) for large N,M. One finds: 
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s- — f :-!].,] )+0(a 2 s ) 
( 21 ) 

We note that an alternative, but equivalent, form of 
the resummed result is [TH HE] 
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This expression can be obtained from ( 15 ) by first writing 
the integrand as 



e sq 2 - 1 = (c 4 



1 
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1 + c M ^ - 1. (23) 



Including the integrals over £ and k± and substitu ting 
k\ = (£Q 2 , the first term on the right-hand-side of (23) 
yields 



c 



- M <-1) A q {aMQ 2 )) 



(24) 

Next, we deal with the term containing (e N * — l) 
in (23 1. Combining with the logarithm In A in (151 we 
find. 



, up to corrections suppressed as 1/N: 
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Adding the expressions in Eqs. ( 24 1 , ( 25 ) , ( 26 ) , and chang- 



ing £ — ► I — £, £ — > I — £, we recover (22). Our expres- 



sion ( 20 ) is however simpler and makes the close connec- 



tion to the resummed expression for the total Drell-Yan 
cross section more transparent. 



C. Expansion to NLL 

At small argument, the Bessel function K behaves as 
K {x) = - In {xe' E /2) + 0(x 2 Inx) . 



Therefore, logarithmic behavior of the integral in (20) in 
TV and M occurs only when fcj^ is bounded from below, 
and the K$ function effectively acts as step function. To 
NLL, one finds the condition k ± > Q/VNM [IT]: 
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(27) 



The explicit NLL expansion of the right-hand side is 
straightforward after inserting the standard expression 
for the running strong coupling and, anyway, the result 
can be directly obtained from the well-known one [26] 



for Drell-Yan, setting N — > V NM there. In the expo- 



nent in (27) we obtain 
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where 



X NM = b Q a s (fi 2 ) (log TV + log M) 



2Trb X 
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with 
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(30) 



The functions h q x \ h q 2 ^ collect all leading-logarithmic 
and NLL terms in the exponent, which are of the form 
oc s In™ iVln m M with n + m = k + 1 and n + m = k, 
respectively. Note that we have kept the factorization 
and renormalization scales arbitrary in the above expres- 
sions. The standard Drell-Yan result is recovered by set- 
ting AjvAf/2 A D y, where A D y = b a s (fi 2 ) log N. 



D. Inverse Mellin transforms 

As the exponentiation of soft-gluon corrections is 
achieved in Mellin moment space, the hadronic structure 
function is obtained by taking the inverse Mellin trans- 
forms of Eq. ([8|: 



dN 

2~TTi' 



-N 



dM 
2~7ri 



z- IU T?(N,M,Q'), 
(31) 

where Cn and Cm denote integration contours in the com- 
plex plane, one for each Mellin inverse. When performing 
an inverse Mellin transform, the contour usually has to 
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be chosen in such a way that all singularities of the in- 
tegrand lie to its left. However, as can be seen from 
Eq. (28), (29), the resummed cross section has a Landau 
singularity at Xnm = 1 or 



are exponentially suppressed by the factors x 



— N „—M 
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3 l/(a s bo)-27E 



(32) 



as a result of the divergence of the running coupling a s 
in Eq. (20) for k± —> Aqcd- For the Mellin inversions, 



we adopt the minimal prescription developed in Ref. [26] 
to deal with the Landau pole. For this prescription the 
contours are chosen to lie to the left of the Landau singu- 
larity. In order to achieve this, we first choose the contour 
for the iV-integration in the upper complex half plane in 
the standard way (see Fig. [TJa)) as 



N = cn + zn e 



i<t>* 



(33) 



where cjv is a postive real constant, <j>N ~ 3ir/4 and 
zjv G [0, oo]. For the lower branch of the contour, one 
simply uses the complex conjugate of TV. The contour for 
the M-integration is parameterized in a similar fashion, 
with a constant cm, an integration parameter zm and 
an angle </>m that we will address shortly. The fact that 
both contours are tilted into the half-plane with nega- 
tive real part improves the numerical convergence of the 
integration, since contributions with negative real part 



Eq. (31) 



In accordance with the minimal prescription [26] , the 
parameter cn is chosen to be smaller than Nq = Lq/cm- 
As N moves along its contour from a point with large 
negative imaginary part to a point with large positive 
one, the Landau pole given by Eq. p2| ) describes a tra- 
jectory shown by the solid line in Fig. |TJb) . The angle 
<$M is now chosen in such a way that the M-integration 
contour (shown by the dashed line) avoids this trajec- 
tory. The larger the imaginary part of N, the larger 
the angle 4>m needs to become. Of course, ultimately as 
3(iV) — > ±oo, the Landau pole moves to the origin in 
the M-plane, and the contour falls onto the real M-axis. 
However, as described above, the contributions from such 
large values of N are extremely suppressed. We note that 
a similar approach to the choice of contours was discussed 
in Ref. |27) . where combined inverse Mellin and Fourier 
transforms were considered. An alternative approach is 
to expand the resummed formula to high perturbative or- 
ders. At each finite order, the Landau pole is not present 
and, therefore, standard Mellin contours can be chosen. 

We match the resummed cross section to the NLO one 
by subtracting the 0(a s ) expansion of the resummed ex- 
pression and adding the full NLO cross section. This 
"matched" cross section consequently not only resums 
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the large threshold logarithms to all orders, but also con- 
tains the full NLO results for the q —> q, q — > g and g — > q 
channels. We will occasionally also consider a resummed 
cross section that has not been matched to the NLO one. 
We will refer to such a cross section as "unmatched" . 



+ Bq (a s ((l-OQ 2 )) 



(39) 



The NLL expansions of these functions can be performed 
in the same way as described in the previous subsection. 
One obtains: 



E. Resummation for inclusive DIS 

In order to obtain resummed predictions for the SIDIS 
hadron multiplicities defined in ([3]), we also need the re- 
summation for the inclusive cross section. In this case, 
there is only one variable x, and standard Mellin- moment 
resummation techniques |28[ 129] may be applied. The 
inclusive structure functions Ft and Te introduced in 
Eq. Q can be written as 

Q 2 ) = £ £ f / (f , „ 2 ) C) (x, J, aV) 

(34) 

We refer the reader to Ref. [21] for the NLO expressions 
for the coefficient functions Cj. Resummation may be 
performed by again introducing Mellin moments in x. 
The resummed DIS coefficient function for the structure 
function Tt reads to NLL in moment space: 



CT' reS = e 2 i? g U(Q 2 ), ., 



exp 



log A? = \nNh^(\) + h^ (x,%,9l), 

V A* A*f / 



(40) 



with the functions h^\h^ given in (29). Furthermore, 



with 



lnjf = lniV/W(A) + /f)(A,^) ! (41) 

8 

/, (1) (A) = Q)-fc«(A), 
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+ ^-ln(l-A). 



III. RESUMMATION FOR e+e" 



(42) 
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^ A q (a s (kl)) + -B q (a s ((l-0Q 2 )) 



(35) 



where the function A q (a s ) is as in ( 16 ). The perturbative 
function B q (a s ) is given by 



with 



B q (a s ) = ^B q ^ +0(a 2 s ), 



Finally, the hard-scattering coefficient reads 



H q ( a s (Q 2 " Q 



(36) 



(37) 
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Hadron multiplicities in e + e — > hX are defined by 

1 d 2 C7 h 



cr tot dx Ed cos 9 



(43) 



where d 2 cr h / dx Ed cos 9 is the differential cross section for 
the production of the hadron h at angle 9 relative to the 
initial positron. Furthermore, 



2P h ■ q 

Q 2 ' 



(44) 



where Ph and q are the momenta of the produced hadron 
and the intermediate virtual photon respectively and 
Q 2 = q 2 . <7 tot is the total cross section for e + e _ — > 
hadrons. To first order in a, it reads: 
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Af x Jf , where 



The exponential in ( 35 ) can evidently be written as 
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logjf ee 
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ds^- r / -vf A,(a fl (fci)) 
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where Ac = 3 is the number of colors. 

As in the case of SIDIS, one can write the cross section 
d 2 a h /dxEd cos 9 in terms of structure functions [2U1 [2Tj : 
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where 
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-Q 2 , 2, 
z, -j ,a a (/x ) 



7r u SIDIS Multiplicity 



(47) 

with the fragmentation functions introduced in sub- 
section |II A| To lowest order, only the partonic channel 
e + e~ —¥ qX contributes, for which 



e q S(l - z), 



C^\z) = 0. 



(48) 



Again we refer the reader to the previous literature [T^l 
|2T| [30] for the NLO expressions for the coefficient func- 
tions. After taking Mellin moments in Xe, the resummed 
result for the corresponding hard-scattering function 
turns out to be identical to that in ( 35 1 (T^HH] , except 
for a change — 7r 2 /6 — > 5ir 2 /6 in the coefficient H q in (38 1. 



IV. PHENOMENOLOGICAL RESULTS 

We now investigate the numerical size of the thresh- 
old resummation effects for the two semi-inclusive 
hadron production processes discussed above, SIDIS and 
e + e~ — > hX. We focus entirely on pion production in this 
work, for which the theory is expected to be under best 
control. We also consider a proton target throughout this 
work. For the parton distribution functions we use the 
NLO "Martin-Stirling-Thorne-Watt" (MSTW 2008) set 
of [31], whereas we choose the NLO "de Florian-Sassot- 
Stratmann" (DSS) [T] pion fragmentation functions. In 
this set, fragmentation functions for charged pions -n^ 
are separately available. We note that the parton dis- 
tributions and fragmentation functions are provided in x 
(or z) space, whereas according to Eq. Q we need their 
Mellin moments. To obtain the latter, we first fit suitable 
functions of the form Ax a (l — x)^ times a polynomial in 
x to the distributions. It is then straightforward to take 
Mellin moments of the fitted functions analytically and 
use them in the numerical code. We have checked that 
the accuracy of the fit is overall very good. 



A. Results for SIDIS 

We start by examining the overall effects of threshold 
resummation for SIDIS, using the kinematics relevant for 
the COMPASS SIDIS measurements [5] as an example. 
COMPASS uses a muon beam of energy 160 GeV inci- 
dent on a proton fixed target. The resulting center-of- 
mass energy is yfs «17.4 GeV. The kinematic cuts em- 
ployed by COMPASS are 0.041 < x < 0.7, O.K y < 0.9, 
Q 2 > 1 GeV 2 and W 2 = Q 2 (l - x)/x + to 2 > 49 GeV 2 , 
where m p is the proton mass. We choose the renormal- 
ization and factorization scales as Q and consider the 
SIDIS multiplicity for neutral pions 7r° as a function of 
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FIG. 2: SIDIS multiplicity for 7r°. Kinematics are as for the 
COMPASS measurements [8J (see text) with, however, an addi- 
tional cut x > 0.4 to enhance the contribution from the threshold 
regime. The lower two lines which are almost indistinguishable 
show the NLO result and the first-order expansion of the re- 
summed result (without matching). The top line shows the full 
NLL resummed result; the lines in between display various fixed- 
order expansions of the latter. 



z, but integrating numerator and denominator of Eq. |3]) 
over x and y. We observe that the range of x probed by 
the full COMPASS data sample extends down to fairly 
low values, where one could be quite far from the thresh- 
old regime. We therefore first consider a lower cut of 
x > 0.4. Figure [2] shows our results for the NLO and the 
resummed multiplicities, along with those for expansions 
of the resummed cross section to various finite orders in 
a s . In each case, the denominator of the multiplicity, the 
inclusive-DIS cross section, has been treated in the same 
fashion as the SIDIS one in the numerator. With the 
exception of the first-order expansion, all beyond-NLO 
results have been matched to the NLO one (separately 
in the numerator and the denominator) as described at 
the end of Sec lIIDI We first of all note that the first-order 
expansion of the resummed cross section agrees very well 
with the full NLO one, which demonstrates that for the 
chosen kinematics the threshold regime strongly domi- 
nates the NLO cross section. The full resummed result 
shows a marked increase over the NLO one, in particu- 
lar at high z, and the higher-order expansions converge 
nicely to the resummed result. Clearly, the a 2 , and cv^ 
contributions generated by resummation are still signifi- 
cant. 

The results can also be studied as ratios 
(Th' - NLO)/NLO, where Th' denotes any of the 
higher-order SIDIS multiplicities generated by resum- 
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z 



FIG. 3: Ratios (TV - NLO)/NLO, where Th' corresponds to 
the SIDIS multiplicity at higher orders as generated by resum- 
mation. The first-order expansion of the resummed cross section 
has not been matched since it would otherwise be identical to 
NLO. 

mation. Figure [3] shows these ratios for the unmatched 
first-order expansion, the matched higher-order expan- 
sions and the full resummed result. The good agreement 
of NLO and the first-order expansion is evident, as are 
the large resummation effects at high z. 

We now extend the x-range to the full region 0.041 < 
x < 0.7 covered by COMPASS. Figs. [4] and [5] show the 
corresponding results, where all lines directly correspond 
to the ones shown in Figs. [2] and [3] for the case x > 0.4. 
One can see that the resummation effects arc generally 
smaller now, even though they remain significant at high 
z. As expected, the agreement between NLO and the 
first-order expansion of the resummed cross section is 
worse now, but it typically remains at the 10% level or 
better. The second-order expansion captures most of the 
full resummation effects; the yet higher orders converge 
somewhat more slowly now to the resummed result. Fig- 
ure[5]shows the corresponding ratios (Th' — NLO)/NLO, 
where again Th' denotes any of the SIDIS multiplicities 
computed at higher orders. 

Preliminary precise data from COMPASS for charged- 
pion and kaon multiplicities are available for a wide range 
of kinematics [5j. The full data set, which evidently has 
the best statistics, covers the range used above for Figs. [4] 
and [5j Figures [6] and [7] show comparisons of our NLO 
and NLL resummed calculations for charged pions to the 
COMPASS data. As one can see, resummation leads to 
a moderate, but significant, enhancement of the multi- 
plicities. It is interesting to note that such an enhance- 
ment is in fact preferred by the ir~ data. However, we 
do not assign much importance to this observation. The 
fragmentation functions are presently still not very well 



7T° SIDIS Multiplicity 




0.2 0.4 0.6 0.8 
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FIG. 4: Same as Fig. [2] but for x > 0.041. 
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FIG. 5: Same as Fig. [3] but for x > 0.041. 



determined and a different set (or a new fit) might well 
describe the data also at NLO. Our main point is that 
inclusion of resummation effects in an analysis of frag- 
mentation functions could make a significant difference 
for the extracted functions. 

Figures [§] and [9] show similar comparisons to the HER- 
MES preliminary data [TJ. For this data set, the center- 
of-mass energy is y/s «17.4 GeV. The kinematic cuts em- 
ployed by HERMES are 0.023 < x < 0.6, O.K y < 0.85, 
Q 2 > 1 GeV 2 and W 2 = 10 GeV 2 . The results are 
qualitatively similar to those shown for COMPASS kine- 
matics. 
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FIG. 6: NLO and NLL resummed SIDIS multiplicities for tt+. 
The results are compared to the preliminary COMPASS data [8.. 
The uncertainties of the data are smaller than the symbol size 
used in the plot. 



FIG. 8: Same as Fig. [6] but for HERMES kinematics. The 
preliminary data are from [7J. The uncertainties of the data are 
smaller than the symbol size used in the plot. 
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FIG. 7: Same as Fig. [6j but for tt~ 
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FIG. 9: Same as Fig. [8j but for tt~. 



B. Results for single-inclusive e + e annihilation 



Figure 10 presents our results for the n~ multiplic- 
ity in e + e~ annihilation at y/s = 10.52 GeV and for 
— 1 < cos 9 < 1, as appropriate for comparison to the 
forthcoming BELLE data [5]. The tt + multiplicity is 
identical thanks to charge conjugation symmetry. We 
have chosen the factorization and renormalization scales 
as -J~s. As for SIDIS. we show NLO and NLL resummed 



results, along with various fixed-order expansions of the 
resummed multiplicity. One first of all observes the 
excellent agreement between the NLO result and the 
first-order expansion of the resummed one. This clearly 
demonstrates that the threshold regime strongly dom- 
inates for the BELLE kinematics. We can therefore be 
confident that also the resummed result reliably captures 
the important higher-order terms. Resummation leads to 
a significant enhancement of the 7r~ multiplicity. This en- 
hancement becomes particularly strong at high xe, but 



11 



77 Multiplicity in e + e 



X 1 
T3 10 1 



10 u 



^ 10" 



10" 



10" 



10" 



T 



T 



1 1 1 I 1 1 1 

VS= 10.52 GeV 



Resummed 
Expansion 5th order 
Expansion 3rd order 
Expansion 2nd order 
Expansion 1st order 
NLO 



0.2 



0.4 



0.6 



0.8 



FIG. 10: tt~ multiplicity in electron-positron annihilation at 
^/s = 10.52 GeV. The two lowest lines show the NLO result 
and the (unmatched) first-order expansion of the resummed one, 
which are practically indistinguishable. The other lines show 
matched higher-order expansions of the resummed multiplicity, 
and the resummed result itself (solid line). 



is present also at moderate values. Figure [TT] shows the 
corresponding ratios (Th' - NLO)/NLO, where Th' de- 
notes any ol the higher-order multiplicities generated by 
resummation and shown in Fig. |10| Since the precision 
of the preliminary BELLE data is typically much better 
than 10%, it will be important to include the enhance- 
ments we hnd in future global analyses of fragmentation 
functions, similar to what has been done in the past in 
Ref. [3J. 



V. CONCLUSIONS 

We have derived threshold-resummed expressions for 
the coefficient functions for single-inclusive hadron pro- 
duction in semi- inclusive lepton scattering and e + e~ an- 
nihilation. We have presented phenomcnological results 
for pion multiplicities for these processes in the kinematic 
regimes presently accessed by the COMPASS, HERMES 
and BELLE experiments. We have found that resumma- 
tion leads to modest but significant enhancements of the 
multiplicities. The recent preliminary SIDIS data turn 
out to be overall better described when resummation ef- 
fects are included in the calculation, at least for the DSS 
set of fragmentation functions that we have used. How- 
ever, we do not ascribe much significance to this point as 
the fragmentation functions are still rather poorly con- 
strained so that a new NLO fit to the new data would 
likely also work well [2], Our main point is that, given the 
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FIG. 11: Ratios (Th' - NLO)/NLO corresponding to the vari- 
ous curves shown in Fig. [To] 



good accuracy of the new preliminary SIDIS and BELLE 
data, it will be crucial to include resummation effects for 
both processes in the next generation of analyses of frag- 
mentation functions. In fact, because of the enhance- 
ments from threshold resummation, the extracted frag- 
mentation functions would be expected to be softer or 
smaller at high z than functions extracted purely on the 
basis of an NLO framework. This may well have impor- 
tant ramifications for the QCD predictions obtained for 
other processes sensitive to fragmentation functions. For 
instance, it has recently been observed [3] [32] that the 
ALICE data [33J for neutral-pion production at 7 TeV are 
well below the theoretical NLO expectations. One may 
speculate if this is due in part to fragmentation functions 
that are too large at high z. 

We note that at very high z or xe, besides the log- 
arithmic perturbative corrections also nonperturbative 
power corrections will ultimately become relevant and 
will need to be analyzed theoretically. Resummation of- 
fers ways to address these contributions (see, for exam- 
ple [33] , and references therein) . Based on the ideas 
presented there, we do not think that power corrections 
play an overwhelming role in the presently accessible 
kinematic regime in SIDIS. As we saw in Eq. (27), the 



logarithmic contributions to SIDIS come from the region 
kj_ > Q/VWm. In ir-z-space this scale roughly corre- 
sponds to QyJ (1 — z)(l — x). We have checked that for 
COMPASS kinematics even at z = 0.9 the average value 
of this scale is significantly larger than 1 GeV, implying 
that perturbation theory should still provide a reliable 
answer here. This issue obviously deserves a more de- 
tailed investigation in the future, in particular also for 
the case of e + e" annihilation. 
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We stress that our study may be extended in several 
ways. First, as mentioned in the Introduction, resum- 
mation for e + eT — > hX could be carried out at next- 
to-next-to-leading logarithm and even beyond, thanks 
to [T31 H3] . The high precision of the BELLE data may 
well warrant a future study along these lines. Resum- 
mation for SIDIS could probably also be extended to 
next-to-next-to-leading logarithmic accuracy with some 
moderate further developments. This may become a very 
worthwhile task in the future when high-precision SIDIS 
data will become available from measurements at the Jef- 
ferson Laboratory after the CEBAF upgrade to 12 GeV 
beam energy |35j . Finally, one may readily adapt the 
rcsummation framework to the case of polarized SIDIS, 
which serves as an important probe of the spin structure 
of the nucleon. Present polarized-SIDIS measurements 
have the same kinematics as those considered in this pa- 
per. 
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Appendix 

In this Appendix we collect the NLO expressions for 
the spin-averaged partonic SIDIS cross sections. At 
NLO, we have to consider the processes j*q — > qX, 
7*q — > gX, and 7*3 —> qX. For their contributions to the 
structure function Tp we have in the MS scheme 



C^(x,z) = e 2 q C t 



(l-x)5(l-z) 



+6(1 -x) 
+6(1 - z) 
+ 



P qq (z) In % + L x (z) + L 2 (z) + (l-z) 



P qq (x) ln-%- +L 1 (x) -L 2 (x) + (l-x) 



1 + z 



1 + x 



(l-x) + (l-z) + (l-x) + (l-z) A 
+2(1 + xz) 



(49) 



CT^(x,z) = elC B 



P gq (z)[5(l~x)\n[^z(l-z) 



Q 2 



(1-x) 



1 \ 1 + x 
+ z6(l -x) + 2(l+x- xz) — 



(50) 



6(1 -z) 



P qg (x)\n[Q T l 



Hp x 



+2x(l-x) +P qg (x) 



(l-z) + z 



, (51) 



where e q is the quark's fractional charge, Cp = 4/3, Tr 
1/2, 





1 


+e 3 


(1 






_ 1- 


-a-o 2 




a 




= e 


+ (i-o 2 , 


Li(0 - 


= (U 




L 2 (0 -- 


_ 1 + 
1 - 





(52) 



and the "+" - distributions are defined as follows: 







«(/(0-/(i))ff(0. 



JO 



f(x,z) 



1 ,1 

dx I dz- 

Jo I 1 -^J+l 1 -*)+ 



1 , z f(x,z)-f(l,z)-f(x,l) + f(l,l) 

(53) 



= / dx dz 

Jo Jo 1-5(1-*) 



Note that we have given expressions ( 49 ]- ( 51 1 for an ar- 
bitrary factorization scale fip, keeping however the scales 
the same for the initial and the final state. For the lon- 
gitudinal structure function F^: 

C^ q {1) (x,z) = 4e 2 q C F xz, 

C^ (1 \x,z) = 4e 2 q C F x(l-z), 

C%W{x,z) = 8e 2 q T R x(l-x). (54) 
In Mellin-moment space, the NLO results become [25] 

1 2 1 

AP + (M+ l) 2 + JP 



C T J l \N,M) = e 2 q C t 



-8- 



(1 + M + A) 2 + 1 



■ 3S 2 (M) - S 2 (N) 



M(M + 1)N(N + 1) 
+ [Si(M) + Si (AO] \ Si(M) + 5i(JV) 



1 1 

' M(M + 1) ~ N(N + 1) 
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—r r + 3 -2SAN) 

N(N + 1) 2 v ' 



In 



where 



1 



M(M + 1) 



25i(M) 



In 



91 



(55) 



(58) 



Note that at large N we have 



(7£W(iV,M) 



e 2 CV 



2M 



2-2M-9M 2 + M 3 ~M 4 + M 5 Sl (jV) = lniV + C^l/iV), S 2 (JV) = — + <D(l/N), (59) 
M 2 (M — 1) 2 (M + l) 2 6 

where N — Ne lB . Furthermore, for the longitudinal 



2 - M + M 2 



structure function, 



N(M + 1)(M - 1) M(M + 1)(M - 1)(JV + 1) 
2 4- M 4- M 2 



2 + M + M 2 /Q 2 
M(M + 1)(M-1) n \n F 

2 + N + N 2 



C%W(N,M) = e 2 q C F 
C^\N,M) = e 2 q C F 



(56) 



C^\N,M) = e 2 Q Tj 



q r (Af+l)(JV+l)' 
4 

M(M + 1)(N + 1) 
8 



&?M(N,M) = e 2 T 



q M(N + l)(N + 2) 



1 Here we have corrected a mistake in Cqg(N,M) in 



(60) 

, (61) 
(62) 



Ar(JV + l)(iV + 2) \ M- 1 M Ref. [25]. 



-5i(M)-5i(JV) + ln 
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N 2 
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